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Abstract 

The  objective  of  this  research  is  to  estimate  the  micro  structure  dependence  of  the  extreme  value  probabilities  of  fatigue  crack 
formation  at  the  scale  of  the  grains  in  polycrystalline  and  polyphase  microstructures  to  facilitate  preliminary  parametric  design 
exploration  and  property  assessment.  A  simulation-based  methodology  is  introduced  for  computing  correlation  functions  of 
micro  structure  attributes  marked  by  the  extreme  value  micro  structure- scale  fatigue  indicator  parameters.  Multiple  statistically 
representative  micro  structure  volumes  simulated  using  a  crystal  plasticity  constitutive  model  coupled  with  the  FEM  are  used  to 
characterize  these  extreme  value  marked  correlation  functions.  By  comparing  these  marked  correlation  functions  to  traditional 
correlation  functions  sampled  from  the  bulk  material,  the  interacting  micro  structure  attributes  important  to  the  extreme  value 
response  are  identified.  This  methodology  is  applied  to  identify  important  crystallographic  relationships  relative  to  fatigue  crack 
formation  in  the  polycrystalline  P/M  Ni-base  superalloy  INI 00  and  duplex  Ti-6A1-4V.  It  was  demonstrated  that  certain  high 
contrast  orientation  relationships  in  both  IN  100  and  Ti-6A1-4V  (where  grains  oriented  for  hard  slip  are  near  grains  oriented  for 
easier  slip)  can  significantly  increase  the  driving  forces  for  fatigue  crack  formation. 

Keywords:  metal  fatigue;  extreme  value  statisties;  marked  eorrelation  fiinetions;  mierostrueture-sensitive  design;  Ni-base  superalloys;  Ti-6A1-4V. 


1.  Introduction 

Component  service  life  under  cyclic  loads  (i.e.  fatigue)  is  primarily  a  function  of  the  extreme  value  probabilities 
of  life  limiting  micro  structure  attributes  or  combinations  of  micro  structure  attributes.  The  life  limiting 
microstructure  attributes  are  those  that  increase  the  driving  forces  for  fatigue  crack  formation  the  most  within  a 
given  volume  of  material.  Moreover,  multiple  interacting  micro  structure  attributes  that  exist  together  with  varying 
probabilities  can  further  enhance  local  driving  forces  for  fatigue  crack  formation.  Traditionally,  the  extreme  value 
microstructure  dependence  of  components  subjected  to  cyclic  loads  has  only  been  captured  empirically  through 
extensive  experimentation.  In  the  last  few  decades,  much  work  has  been  done  to  better  understand  the 
microstructure  dependent  mechanisms  of  fatigue  crack  formation  in  metals  [1,2].  However,  to  predict  variability  in 
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the  response  of  damage  dependent  processes  like  fatigue  crack  formation,  a  framework  is  needed  that  accounts  for 
the  extreme  value  life  limiting  attributes  most  important  to  the  processes  of  fatigue  crack  formation. 

Current  probabilistic  approaches  that  account  for  the  extreme  value  statistics  of  life  limiting  micro  structure 
attributes  relative  to  fatigue  crack  formation  consider  only  distributions  of  a  single  micro  structure  attribute  and  do 
not  account  for  the  effect  that  multiple  interacting  micro  structure  attributes  has  on  the  extreme  value  response 
distributions.  For  example,  Murakami  et  al  [3,  4]  use  classical  extreme  value  Gumbel  statistics  to  estimate  the  size 
of  the  largest  inclusion  based  on  specimen  volume.  In  this  case,  the  authors  use  extreme  value  statistics  of  inclusion 
size  to  estimate  the  fatigue  strength  in  clean  steels  based  on  the  assumption  that  the  life  limiting  fatigue  crack  tends 
to  form  at  the  largest  inclusion  within  a  given  volume,.  Similarly,  Atkinson  and  Shi  [5]  use  a  generalized  Pereto 
distribution  to  estimate  the  size  of  the  largest  inclusion  for  a  given  volume.  The  difference  being  that  the  generalized 
Pereto  distribution  incorporates  a  limit  on  the  maximum  inclusion  size;  whereas,  the  extreme  value  Gumbel 
distribution  predicts  an  ever  increasing  inclusion  size  for  increasing  volumes.  These  types  of  models  based  on  a 
single  microstructure  attribute  like  inclusion  size  do  not  account  for  the  interactions  between  the  inclusion  and 
surrounding  matrix,  for  example.  Although  first  order  approaches  based  on  single  micro  structure  attributes  may  be 
sufficient  for  some  material  systems  or  applications,  the  process  of  fatigue  crack  formation  in  many  advanced 
engineering  alloys  can  be  very  complex.  As  such,  to  understand  the  processes  of  fatigue  crack  formation  in  complex 
material  systems,  a  better  probabilistic  description  of  the  interactions  between  multiple  micro  structure  attributes  is 
required. 

The  objective  of  this  research  is  to  apply  a  newly  developed  extreme  value  micro  structure  sensitive  probabilistic 
framework  to  quantify  the  influence  of  polycrystalline  micro  structure  on  fatigue  crack  formation  in  the 
polycrystalline  P/M  Ni-base  superalloy  INI 00  and  in  duplex  Ti-6A1-4V  under  HCF  conditions.  Specifically,  this 
new  probabilistic  framework  uses  extreme  value  marked  correlation  functions  for  micro  structure  attributes  to 
quantify  the  interacting  micro  structure  attributes  that  are  most  statistically  significant  relative  to  the  measured  or 
simulated  extreme  value  distributions  of  response  [6,  7].  Correlation  functions  for  the  micro  structure  attributes 
important  to  the  operant  mechanisms  of  fatigue  crack  formation  are  selectively  sampled  in  regions  marked  by 
certain  extreme  value  response  parameters  {e.g.  cyclic  plastic  strain  range).  These  extreme  value  marked 
correlations  can  then  be  compared  to  the  same  correlation  functions  sampled  over  the  entire  micro  structure 
ensemble  (i.e.,  without  regard  to  response)  to  determine  the  particular  correlated  micro  structure  attributes  that  are 
most  significant.  In  general,  this  methodology  can  be  applied  to  quantify  important  correlated  micro  structure 
attributes  relative  to  their  influence  on  the  extreme  value  distributions  of  response  for  many  different  damage 
processes  (e.g.  fatigue,  creep,  fracture).  In  this  work,  multiple  material  volumes  for  both  INI 00  and  Ti-6A1-4V  are 
simulated  via  the  FEM  with  crystal  plasticity  constitutive  relations  to  explore  how  local  crystallography  influences 
the  distributions  of  the  estimated  driving  forces  for  fatigue  crack  formation  as  estimated  by  certain  cyclic  plasticity 
based  fatigue  indicator  parameters  (FIP). 

1.1.  Fatigue  crack  formation  in  poly  crystalline  P/M  Ni-base  superalloy  IN  100 

Typically,  fatigue  damage  formation  in  polycrystalline  superalloys  has  been  linked  to  the  existence  of  large  pores 
or  nonmetallic  inclusions  introduced  during  processing.  Fatigue  cracks  form  at  inclusions/pores  under  cyclic  loading 
because  of  local  stress  concentrations  [8-10].  In  some  cases,  however,  cleaner  polycrystalline  Ni-base  superalloys 
are  being  developed  that  exhibit  fatigue  crack  formation  along  crystallographic  planes  in  samples  or  components 
absent  of  inclusions/pores.  For  example,  Jha  et  al.  [11]  noted  that  subsurface  fatigue  crack  formation  occurs  in 
individual  grains  absent  of  any  voids/inclusions  in  Rene  88DT  in  the  HCF  regime  at  593  °C.  Here  the 
crystallographic  facets  associated  with  fatigue  crack  formation  were  observed  to  be  much  larger  than  the  average 
grain  size,  suggesting  that  cracks  tend  to  form  in  larger  grains.  Similar  results  were  reported  by  Miao  et  al.  [12]  in 
the  same  superalloy  in  the  VHCF  regime.  In  particular,  Miao  et  al.  noted  that  the  faceted  grains  near  the  fatigue 
crack  formation  site  were  associated  with  large  grains,  higher  Schmid  factors  for  {1  1  1}<1  1  0>  slip  and  twin 
boundaries  of  a  particular  character. 

Different  slip  modes  appear  to  be  more  important  depending  on  the  service  temperature.  The  y  austenitic  phase  in 
Ni-base  superalloys  has  a  fee  crystalline  lattice,  which  has  12  octahedral  {11  !}<!  1  0>  slip  systems.  Moreover  at 
certain  elevated  temperatures  an  additional  6  cube  {1  00}<1  1  0>  slip  systems  have  been  observed  in  the  y  phase  in 
many  Ni-base  superalloys  in  the  form  of  macroscopic  slip  traces  [13].  Although  limited  cube  slip  is  possible  in  the 
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LI2  ordered  y'  phase,  it  is  not  commonly  considered  as  likely  to  occur  in  the  austenitic  y  phase.  Based  on 
observations  via  transmission  electron  microscopy  (TEM),  Bettge  and  Osterle  [13],  postulated  that  cube  slip 
observed  along  {10  0}  planes  is  actually  due  to  the  “zig-zag”  cross  slip  of  screw  dislocations  on  {1  1  1}  planes 
when  blocked  at  the  y/  y’  interfaces.  As  these  dislocation  migrate  within  the  channels  between  y'  precipitates,  they 
create  macroscopic  {10  0}  slip  traces.  In  the  polycrystalline  P/M  Ni-base  superalloy  INI 00,  fatigue  cracking  has 
been  observed  to  propagate  along  {111}  planes  at  lower  homologous  temperatures  and  along  {10  0}  planes  at 
moderate  temperatures  (i.e.,  between  500°C  to  800°C)  [14]. 

7.2.  Fatigue  crack  formation  in  a+f  Ti  alloys 

Fatigue  crack  formation  in  a+y?  Ti  alloys  is  primarily  associated  with  the  development  of  crystallographic  facets 
at  the  grain  scale.  Bridier  et  al  [15,  16]  investigated  fatigue  crack  formation  in  Ti-6A1-4V  in  both  LCF  and  HCF, 
respectively,  and  found  that  cracks  primarily  form  on  basal  slip  planes  and  less  frequently  on  prismatic  slip  planes 
favorably  oriented  for  slip.  Gilbert  and  Piehler  [17]  observed  internal  fatigue  crack  formation  on  pyramidal  planes  in 
the  primary  a  grains  in  Ti-6A1-4V  oriented  for  hard  basal  slip  (i.e.,  with  c  axis  nearly  parallel  to  the  loading 
direction).  Clusters  of  grains  oriented  for  slip  on  pyramidal  <a+c>  planes  and  grains  oriented  for  basal  slip  were 
observed  by  the  authors  to  enhance  fatigue  crack  formation.  Szczepanski  et  al  [18]  observed  fatigue  facet  formation 
in  the  VHCF  regime  of  Ti-6246  on  or  near  basal  or  prismatic  slip  planes  favorably  oriented  for  slip.  In  the  same 
material,  Jha  and  Larsen  [19]  noted  that  at  moderate  stress  amplitudes  (i.e.,  close  to  the  macroscopic  yield  stress) 
that  facets  near  the  location  of  initial  fatigue  crack  formation  were  all  inclined  -35-45°  relative  to  the  loading  axis. 
In  one  sample,  the  facet  near  the  location  of  initial  fatigue  crack  formation  appeared  to  be  -35°  relative  to  the 
loading  direction  in  a  grain  oriented  unfavorably  for  basal  slip,  suggesting  the  facet  was  formed  by  a  <c+a>  hard 
slip  mechanism  (i.e.,  pyramidal  slip).  In  another  sample,  the  facet  near  the  location  of  initial  fatigue  crack  formation 
developed  on  an  apparent  basal  slip  plane  favorably  oriented  for  slip.  The  grain  unfavorably  oriented  for  basal  slip 
that  developed  the  facet  was  surrounded  by  other  grains  more  favorably  oriented  for  easy  slip  (i.e.,  prismatic  or 
basal  slip).  Additionally,  the  sample  with  the  facet  that  formed  in  the  grain  oriented  for  hard  basal  slip  failed  at 
-60,000  cycles;  whereas,  in  the  other  sample  with  the  facet  that  formed  in  the  grain  oriented  for  easy  basal  slip 
failed  at  -2,000,000  cycles. 


1.3.  Fatigue  indicator  parameters 


Simulation-based  strategies  for  investigating  the  microstructure-dependence  of  the  driving  forces  for  fatigue 
crack  formation  have  greatly  expanded  in  recent  years.  McDowell  [20]  outlined  the  use  of  plastic  shear  strain-  based 
FIPs  to  estimate  the  driving  forces  for  fatigue  crack  formation  in  metals  due  to  plastic  ratcheting  and  reversed  cyclic 
plasticity.  For  example,  the  maximum  plastic  shear  strain  range  (MPSS)  and  Fatemi-Socie  (FS)  [21]  damage 
parameter,  i.e.. 


P  =- 


yP  ' 

(  max  / 


(1) 


and 


P  = 

^  FS 


r  n  ^ 


(2) 


are  functions  of  the  maximum  plastic  shear  strain  range  The  stress  normal  to  the  plane  of  maximum  plastic 

shear  is  given  by  ay  is  macroscopic  cyclic  yield  stress,  and  K  controls  the  effect  of  the  normal  stress  on  fatigue 
crack  formation.  Calculations  of  the  FS  FIP  for  both  IN  100  and  Ti-6A1-4V  are  carried  out  with  a  value  of  K=\  based 
on  previous  studies  (e.g.  [7]). 

Averaged  over  scales  pertinent  to  the  microstructure  attributes  relevant  the  processes  of  fatigue  crack  formation 
(i.e,  specified  by  F),  these  types  of  FIPs  have  been  used  to  investigate  the  influence  of  several  different 
microstructure  attributes  relative  to  their  potency  for  fatigue  crack  formation.  In  martensitic  gear  steels, 
Prasannavenkatesan  et  al.  [22]  used  these  two  parameters  to  characterize  the  relative  potency  of  different  kinds  of 
inclusions  to  form  grain/inclusion  scale  HCF  cracks.  Similar  FIPs  have  also  been  used  to  characterize  the  influence 
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of  microstmcture  on  fatigue  crack  formation  in  structural  steel  [23,  24]  and  Ni-base  superalloys  [6,  7,  25-27].  These 
efforts  have  demonstrated  the  utility  in  using  these  types  of  parameters  coupled  with  crystal  plasticity  simulations  to 
correlate  microstmcture  scale  slip  with  fatigue  crack  formation  and  early  stages  of  micro stmcturally  small  crack 
growth. 

2.  Methodology 

2.1.  Extreme  value  marked  correlation  functions 

A  new  microstmcture-sensitive  extreme  value  probabilistic  framework  has  been  proposed  that  links  correlated 
microstmcture  attributes  to  extreme  value  response  via  extreme  value  marked  correlation  functions  [6,  7].  To 
constmct  the  extreme  value  marked  correlation  functions,  one  must  first  construct  the  extreme  value  distributions  of 
the  response  parameter  of  interest  {e.g.,  stress,  strain,  plastic  strain).  Based  on  classical  extreme  value  probability 
theory  {cf.  [28]),  the  probability  distribution  of  the  extreme  valued  response  parameter  a  {e.g.,  a=FIP)  is  defined  as 
i^^(a|Q)  given  a  window  Q.  Thus,  F^^{a\Q)  describes  the  probability  that  the  response  parameter  of  value  a  is  the 
extreme  value  for  a  sampled  volume  of  microstructure  Q.  The  extreme  value  marked  radial  correlation  function  as 
^max(a)^  y  is  defined  as  the  probability  that  a  microstmcture  attribute  f  located  coincident  with  the  extreme 
value  response  parameter  a  in  Q  is  located  within  a  distance  of  r  to  r  +  (ir  of  a  second  microstmcture  attribute  f  in 
any  direction.  We  can  similarly  define  an  extreme  value  marked  -point  correlation  function,  extreme  value  marked 
lineal  path  correlation  function,  extreme  value  marked  nearest  neighbor  correlation  function,  etc.  [29].  The  extreme 
value  response  distribution  functions  describe  the  extreme  value  response  of  the  microstmcture  as  represented  by  the 
response  parameter  a.  Coupled  with  extreme  value  response  distribution,  the  extreme  value  marked  correlation 
functions  then  describes  the  probabilities  of  existing  correlated  microstmcture  attributes,  f  and  relative  to  the 
observed  extreme  values  of  a  in  a  microstmcture  window  Q.  Comparing  to  the  correlation 

function  for  the  same  microstmcture  attributes  sampled  from  the  complete  material  ensemble  {e.g.,  bulk)  or 
R(f,fy)dr,  the  correlated  microstmcture  attributes  most  unique  to  the  locations  of  extreme  value  response  can  be 
identified. 

2.2.  Simulated  statistical  volume  elements 

To  characterize  the  microstmcture  attributes  most  important  to  the  slip-dominated  processes  of  fatigue  crack 
formation  in  HCF  in  IN  100  and  Ti-6A1-4V,  here  we  consider  the  crystallographic  attributes  of  the  micro  structure 
{e.g.,  phase,  grain  orientation,  grain  disorientation,  grain  size,  and  grain  shape  distributions).  In  INI 00,  we  recognize 
that  fatigue  crack  formation  is  most  influence  by  the  presence  of  non-metallic  inclusions  and/or  pores.  However,  Ni- 
base  superalloys  have  been  developed  with  much  smaller  populations  of  these  types  of  inclusions/voids  that  fail 
crystallographically  {e.g.,  [30]).  Thus,  here  we  consider  how  this  variant  of  INI  00  would  fail  crystallographically  in 
the  absence  of  non-metallic  inclusions  or  voids.  In  the  case  of  Ti-6A1-4V,  fatigue  crack  formation  is  predominantly 
crystallographic  and  is  very  much  influenced  by  the  microstmcture  attributes  being  considered  here.  The  driving 
forces  for  fatigue  damage  formation  on  the  scale  of  the  grains  is  estimated  by  the  MPSS  and  the  FS  FIP  as  introduce 
previously.  In  each  microstmcture,  these  parameters  are  calculated  over  grain  sized  averaging  volumes. 

To  constmct  the  extreme  value  marked  correlation  functions,  multiple  simulated  statistical  volume  elements 
(SVE)  or  microstmcture  instantiations  (i.e.,  microstmcture  window  Q)  were  constmcted  for  both  INI 00  and  Ti-6A1- 
4V  in  the  form  of  meshed  polycrystalline  micro  structures  with  explicit  grains  and  phases.  Using  the  FEM,  variations 
in  local  crystallographic  response  can  be  simulated  using  crystal  plasticity  constitutive  relations.  Each  SVE  is 
constmcted  based  on  predefined  distributions  of  grain  size,  grain  orientation,  and  grain  misorientation.  Here  we 
define  a  SVE  based  on  the  local  plastic  response.  In  other  words,  each  SVE  is  constructed  to  be  large  enough  such 
that  the  local  plastic  strain  (on  the  scale  of  the  grains)  is  unaffected  by  statistical  variations  in  the  microstmcture  at 
distances  on  the  order  of  the  size  of  the  SVE  [7]. 

The  SVEs  for  IN  100  were  constructed  based  on  a  Voronoi  tessellation  algorithm  that  optimized  the 
microstmcture  to  fit  a  predefine  grain  size  distribution  via  a  simulated  annealing  algorithm  [7].  The  duplex  Ti-6A1- 
4V  SVEs  were  generated  using  an  ellipsoidal  based  microstmcture  generator.  This  method  uses  known  distributions 
of  the  aspect  ratios  of  grain  equivalent  ellipsoids  fit  to  grains  of  experimentally  characterized  micro stmctures.  Using 
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grain  equivalent  ellipsoids  instead  of  other  space  filling  methods  like  Voronoi  tessellation  allows  construction  of 
more  complex  grain  morphologies  such  as  elongated  grains  common  in  rolled  ductile  metals.  For  example,  Brahme 
et  al  [31]  used  ellipsoids  fit  to  actual  grain  distributions  to  reconstruct  A1  micro  structures  with  elongated  grains. 
Groeber  et  al  has  used  similar  methods  to  reconstruct  poly  crystalline  representations  of  IN  100  characterized  from 
three-dimensional  serial  sections  [32,  33].  The  ellipsoid  packing  based  digital  micro  structure  generation  algorithm 
was  implemented  in  MATLAB  [34],  with  logic  as  follows: 

1 .  An  ellipsoid  is  generated  by  randomly  selecting  the  ellipsoid  aspect  ratios  from  a  predefined  distribution 
{e.g.,  beta  distribution). 

2.  The  orientation  of  the  ellipse  is  randomly  selected  based  on  a  predefined  distribution  (i.e.,  orientation 
distribution  function  for  orientation  of  major  axes  of  grain  equivalent  ellipsoids  relative  to  the  specimen 
axis). 

3.  The  ellipsoid  is  randomly  placed  in  the  volume.  If  the  ellipsoid  overlaps  with  any  previously  placed 
ellipsoid,  a  new  random  location  is  selected  repeatedly  until  no  overlap  exists. 

4.  Steps  1  through  4  are  repeated  until  the  jamming  limit  is  reached  and  no  more  non-overlapping 
ellipsoids  can  be  placed  in  the  volume. 

5.  The  ellipsoids  are  then  allowed  to  grow  uniformly  until  all  the  space  in  the  micro  structure  volume  is 
filled. 

A  two-dimensional  section  of  a  generated  micro  structure  using  this  method  with  a  bi-modal  grain  size 
distribution  is  given  in  Figure  1.  The  target  and  generated  bimodal  grain  size  distribution  of  a  SVE  for  Ti-6A1-4V  is 
given  in  Figure  2. 


Figure  1  2-D  section  of  a  microstructure  generated  using  the  ellipsoid 
packing  algorithm  with  a  bi-modal  grain  size  distribution. 


Figure  2  Target  log-normal  and  fit  model  distributions  of  the  grain  size 
for  the  primary  a  grains  and  the  a+p. 


The  crystal  plasticity  constitutive  model  for  IN  100  was  developed  by  Shenoy  et  al  [35,  36].  This  model  does  not 
explicitly  model  the  y’  precipitates,  but  accounts  for  the  influence  of  the  primary,  secondary  and  tertiary  y’ 
precipitates  via  certain  internal  state  variables  incorporated  into  certain  hardening  laws.  This  model  also  accounts  for 
the  average  effect  of  the  grain  size  on  the  constitutive  response  and  the  dependence  of  the  flow  stress  on 
crystallographic  orientation.  Shenoy  et  al  [35,  36]  considers  both  the  standard  12  octahedral  <1  1  0>{1  1  1}  slip 
systems  and  6  cube  slip  systems  {10  0}<1  1  0>.  This  crystal  plasticity  model  for  INIOO  was  calibrated  at  650°C  and 
incorporated  into  a  User  MATerial  Subroutine  (UMAT)  for  use  in  the  commercially  available  FE  package 
ABAQUS  [37]. 

The  crystal  plasticity  model  for  Ti-6A1-4V  was  developed  by  Bridier  et  al  [38]  for  application  to  the  HCF 
loading  regime  for  peak  stresses  at  or  below  the  macroscopic  yield  strength.  The  model  accounts  for  four  different 
families  of  slip  systems  in  the  primary  a  phase:  three  ^1120^(0001)  basal,  three  ^ii20)|loTo}  prismatic,  six 
(i120)|10Ti|  first  order  pyramidal  <a>  and  twelve  (ii23)|l0Tl}  fi^st  order  pyramidal  <a+c>  slip  systems  or  a  total 
of  24  possible  modes  of  slip.  However,  slip  is  dominant  in  the  basal  and  prismatic  slip  systems  due  to  a  relatively 
low  critical  resolved  shear  stress  when  compared  to  the  other  slip  systems.  The  lamellar  a+P  colonies  are 
homogenized  in  this  model.  A  crystallographic  burgers  orientation  relation  (BOR)  is  maintained  between  the 
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secondary  a  and  ^  laths  defined  such  that  (^OOOI)  //{lOl}  (OOOl)  //(lH)  *  possible  slip  systems  in  the 

lamellar  region  include  three  ^ii20\(000l)  basal,  tfiree  ^ll20^|l0T0}  prismatic,^  six  ^ii20^|ioTl|  first-order  pyramidal 
and  twelve  bee  slip  sysiems.  The  bee  slip  systems  are  transformed  into  the  hexagonal  coordinate  system 

via  the  BOR.  Hard  systems  in  these  colonies  are  those  that  intersect  the  a/p  interface.  Soft  deformation  modes  are 
those  on  which  dislocations  glide  parallel  to  the  a/p  interface  or  have  parallel  slip  planes  in  both  the  secondary  a  and 
p  phases.  This  crystal  plasticity  model  was  also  incorporated  as  a  UMAT  in  ABAQUS  [37]. 

In  all  cases,  the  SVEs  were  loaded  in  uniaxial  tension  at  a  quasistatic  strain  rate  of  0.002s'\  Periodic  boundary 
conditions  are  applied  in  all  directions  to  simulate  subsurface  conditions.  The  voxellated  meshes  of  these  generated 
microstructures  consist  of  quadrilateral  elements  with  reduced  integration  (type  C3D8R  in  ABAQUS). 

In  INI 00,  the  number  of  grains  in  each  simulated  SVE  ranged  between  325-375  grains  when  each  SVE  was 
dimensioned  to  0.150  mm  along  each  edge.  The  distribution  of  grain  volumes  normalized  by  the  target  average  grain 
volume  of  8.0x1 0'^mm^  was  fit  to  a  log  normal  distribution  with  a  mean  of  -12  and  standard  deviation  of  0.4.  This 
corresponds  to  a  target  average  cube  root  grain  size  of  ~0. 020mm.  There  are  24  elements  along  each  edge  of  the 
SVE  or  13824  elements  in  all.  For  IN  100,  a  total  of  100  SVEs  were  generated.  The  grain  averaged  FIPs  were 
calculated  after  the  third  simulated  cycle  to  allow  for  shakedown. 

In  Ti-6A1-4V,  the  number  of  grains  in  each  generated  SVE  varied  between  325-375.  The  target  volume  fraction 
of  the  primary  a  phase  was  30%.  The  standard  deviation  of  the  target  log-normal  grain  size  distribution  for  the 
primary  a  phase  was  0.025mm  and  0.010mm,  respectively.  For  the  a+fi  colonies  the  mean  and  standard  deviation 
of  the  target  log-normal  grain  size  distribution  was  0.050mm  and  0.005mm,  respectively.  These  SVEs  were  all 
dimensioned  0.400  mm  along  each  edge.  There  are  26  elements  along  each  edge  of  the  SVE  or  17576  elements  in 
all.  For  Ti-6A1-4V,  a  total  of  100  SVEs  were  generated.  The  grain  averaged  FIPs  were  calculated  after  the  tenth 
cycle. 

3.  Results 

3.1.  Extreme  value  marked  correlations  in  IN  100 


Figure  3  The  extreme  value  grain  averaged  Fatemi-Soeie  (FS)  and 
maximum  plastie  shear  strain  (MPSS)  FIPs  ealeulated  over  100 
simulated  SVEs  for  IN  100  plotted  on  a  Gumbel  probability  seale. 
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Figure  4  The  eross  eorrelations  between  the  Sehmid  faetors  for  both 
the  oetahedral  and  eube  slip  systems  eoineident  with  the  grain  average 
extreme  value  Fatemi-Soeie  FIPs  in  eaeh  simulated  SVE  for  INI 00. 


In  Figure  3  we  plot  the  distributions  of  the  extreme  value  FS  and  MPSS  FIPs  calculated  from  100  simulated  SVEs  of 
INI 00,  each  cycled  three  times  at  0.5%  maximum  applied  strain  under  completely  reversed  loading  (i.e.  E=-l).  This 
plot  is  linearized  for  the  Gumbel  distribution  as  outlined  in  [7].  The  parameters  of  the  Gumbel  distributions  fit  via 
least  squares  regression  to  the  observed  extreme  value  FIPs  are  given  in 
Table  1  where  the  Gumbel  distribution  is  given  as 
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The  characteristic  largest  value  of  the  initial  variable  X  is  given  as  i/„,  and  is  an  inverse  measure  of  dispersion  of 
the  largest  value  of  X.  The  number  of  samples  of  the  initial  variable  X  in  the  set  of  distributions  of  X  from  which 
is  sampled  is  n,  which  for  the  asymptotic  forms  listed  above  is  assumed  to  be  very  large  [39].  The  cross  correlations 
between  the  Schmid  factors  for  both  the  octahedral  and  cube  slip  systems  coincident  with  the  grain  average  extreme 
value  FS  FIPs  are  plotted  in  Figure  4. 


Table  1  Fit  Gumbel  distribution  parameters  for  the  extreme  value  FIP  distributions  for  IN  100. 
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Figure  5  Radial  eorrelation  fimetions  for  the  eomplete  mierostrueture  eompared  with  extreme  value  marked  radial  eorrelation  fiinetions 
deseribing  the  eorrelation  between  the  Sehmid  faetors  m  for  a)  eube  slip  between  0.45  and  0.5  and  oetahedral  slip  between  0.45  and  0.5  and  b) 
eube  slip  between  0.45  and  0.5  and  eube  slip  between  0.45  and  0.5. 

The  extreme  value  marked  correlation  functions  associated  with  the  extreme  value  distributions  of  the  FIPs  are 
constructed  with  the  micro  structure  attribute  of  the  Schmid  factor  (computed  based  on  the  geometry  of  the  slip 
systems  relative  to  the  loading  direction).  By  using  the  Schmid  factor,  the  probabilities  of  correlated  orientations  and 
misorientations  of  slip  systems  are  indirectly  considered.  To  construct  the  extreme  value  marked  correlation 
functions  for  the  Schmid  factor  the  range  of  considered  Schmid  factors  (0.3  to  0.5)  a  tessellated  scale  in  0.05 
increments  is  used.  In  Figure  5,  we  plot  the  complete  radial  distribution  estimated  over  the  entire  ensemble  of 
simulated  SVEs  and  the  extreme  value  marked  distribution  functions  sampled  only  at  the  location  of  the  extreme 
value  FIPs  for  the  100  SVEs  cycled  at  0.5%  strain.  The  particular  correlations  plotted  in  Figure  5  are  for  the  Schmid 
factors  for  a)  cube  slip  between  0.45  and  0.5  and  octahedral  slip  between  0.45  and  0.5  and  b)  cube  slip  between  0.45 
and  0.5  and  cube  slip  between  0.45  and  0.5. 

3.2.  Extreme  value  marked  correlations  in  Ti-6Al-4V 

The  extreme  value  distributions  of  the  FS  and  MPSS  FIPs  are  given  in  Figure  6  for  the  100  SVEs  simulated  for 
Ti-6A1-4V  and  cycled  with  a  positive  mean  stress  (^=0)  to  0.6%  maximum  strain.  The  distributions  of  the  extreme 
value  FIPs  are  plotted  on  a  Gumbel  probability  scale  as  describe  previously.  The  estimated  distribution  parameters 
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for  the  extreme  value  Gumbel  distribution  estimated  via  linear  least  squares  regression  are  given  in  Table  2.  The 
cross  correlations  between  the  Schmid  factors  for  basal,  prismatic,  pyramidal  <a>,  and  pyramidal  <a+c>  slip  and 
the  grain  average  extreme  value  Fatemi-Socie  FIP  estimated  from  each  of  the  100  simulated  SVEs  of  Ti-6A1-4V  are 
plotted  in  Figure  7. 


Extreme  Value  Fatemi-Socie  FIP 

Figure  6  The  extreme  value  grain  averaged  Fatemi-Soeie  (FS)  and 
maximum  plastie  shear  strain  (MPSS)  FIPs  ealeulated  over  100 
simulated  SVEs  for  Ti-6A1-4V  plotted  on  a  Gumbel  probability  seale. 
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Figure  7  The  eross  eorrelations  between  the  Sehmid  faetors  for  basal, 
prismatie,  pyramidal  <a>,  and  pyramidal  <a+e>  slip  and  the  grain 
average  extreme  value  Fatemi-Soeie  FIP  estimated  from  eaeh  of  the 
100  simulated  SVEs  of  Ti-6A1-4V. 


Table  2  Fit  Gumbel  distribution  parameters  for  the  extreme  value  FIP  distributions  for  Ti-6A1-4V. 
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Figure  8  Radial  eorrelation  funetions  for  the  eomplete  mierostrueture  eompared  with  extreme  value  marked  radial  eorrelation  funetions 
deseribing  the  eorrelation  between  Sehmid  faetors  for  a)  basal  slip  between  0.45  and  0.5  for  the  primary  a  phase  and  basal  slip  between  0.45  and 
0.5  for  the  primary  a  phase  and  b)  basal  slip  between  0.45  and  0.5  for  the  primary  a  phase  and  prismatie  slip  between  0.45  and  0.5  for  the 
primary  a  phase. 
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Figure  9  Radial  correlation  functions  for  the  complete  microstructure  compared  with  extreme  value  marked  radial  correlation  ftinctions 
describing  the  correlation  between  Schmid  factors  for  a)  basal  slip  between  0.45  and  0.5  for  the  primary  a  phase  and  (l  1 1  \  1 0}  bee  slip 

(transformed  into  the  hexagonal  coordinate  system  via  the  BOR)  between  0.45  and  0.5  for  the  a+jO  colony  phase  and,  b)  basal  slip  oetween  0.45 
and  0.5  for  the  primary  a  phase  and  pyramidal  (d)  slip  between  0.45  and  0.5  for  the  a+jO  colony  phase. 

As  demonstrated  previously  with  IN  100,  the  extreme  value  marked  correlation  functions  for  Ti-6A1-4V  are 
constructed  using  the  micro  structure  attribute  of  the  Schmid  factor.  As  shown  in  Figure  7,  grains  oriented  for  easy 
basal  slip  were  observed  to  be  predominate  at  locations  of  the  estimated  extreme  value  FIPs.  In  Figure  8,  the  radial 
correlation  functions  for  the  complete  micro  structure  compared  with  extreme  value  marked  radial  correlation 
functions  are  given  describing  the  correlation  between  Schmid  factors  for  a)  basal  slip  between  0.45  and  0.5  for  the 
primary  a  phase  and  basal  slip  between  0.45  and  0.5  for  the  primary  a  phase  and  b)  basal  slip  between  0.45  and  0.5 
for  the  primary  a  phase  and  prismatic  slip  between  0.45  and  0.5  for  the  primary  a  phase.  In  Figure  9,  the  complete 
and  extreme  value  marked  radial  distribution  functions  are  given  for  Schmid  factors  for  a)  basal  slip  between  0.45 
and  0.5  for  the  primary  a  phase  and  1  I^jl  10}  bee  slip  (transformed  into  the  hexagonal  coordinate  system  via  the 
BOR)  between  0.45  and  0.5  for  the  a+jS  colony  phase  and  b)  basal  slip  between  0.45  and  0.5  for  the  primary  a  phase 
and  pyramidal  slip  between  0.45  and  0.5  for  the  a+jS  colony  phase. 


4.  Discussion 


The  extreme  value  FIP  distributions  for  the  FS  and  MPSS  parameters  are  similar  in  character  for  both  INI 00  and 
Ti-6A1-4V.  In  both  material  systems,  there  appears  to  be  less  scatter  in  the  extreme  value  MPSS  FIP  distribution 
than  then  FS  FIP  distribution.  This  can  be  quantified  by  comparing  the  values  of  a„  listed  in  Tables  1  and  2, 
recalling  that  for  the  Gumbel  distribution,  a„  is  an  inverse  measure  of  the  dispersion  of  the  largest  values  of  the 
initial  population.  In  IN  1 00,  the  grains  predicted  to  exhibit  the  extreme  value  response  were  different  between  the 
MPSS  and  FS  FIPs  in  as  many  as  10%  of  the  SVEs  simulated.  There  were  no  grains  with  the  predicted  extreme 
value  response  in  Ti-6A1-4V  with  the  MPSS  FIP  that  were  not  identically  identified  by  the  FS  FIP.  Although  the 
distributions  are  similar  between  the  different  FIPs,  the  FS  FIP  best  captures  the  actual  mechanisms  of  fatigue  crack 
formation  based  on  the  previous  experimental  observations  of  the  local  damage  processes  acting  in  this  class  of 
materials.  As  shown  by  McDowell  and  Berard  [40],  the  FS  parameter  reflects  mixed  mode  growth  of  small  fatigue 
cracks.  For  INI 00  and  Ti-6A1-4V  it  is  expected  to  be  effective  in  reflecting  combined  effects  of  slip  with  elevated 
normal  stress  due  to  intergranular  interactions. 

In  INI  00,  cube  slip  is  predominate  in  the  locations  of  predicted  extreme  value  FIPs  as  evident  from  Figure  4.  In 
Figure  5,  it  is  evident  that  there  is  a  much  higher  probability  of  finding  grains  with  high  Schmid  factors  (i.e.  between 
0.45  and  0.5)  for  cube  slip  near  other  grains  similarly  oriented  or  octahedral  grains  with  similar  high  Schmid  factors 
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at  the  locations  of  the  extreme  value  FS  FIPs  than  in  the  overall  micro  structure.  These  observations  support  the 
results  of  previous  experiments  of  fatigue  in  this  material  system  [41].  It  is  likely  that  these  types  of  orientation 
relationships  with  a  high  contrast  induce  additional  localize  plasticity  as  the  softer  phases  shed  load  onto  the  harder 
phases  during  cyclic  deformation.  Similar  high  contrast  orientation  relationships  were  observed  to  be  important  in 
Ti-6A1-4V. 

In  Ti-6A1-4V,  all  of  the  grains  in  which  the  extreme  value  FIPs  are  identified  were  primary  a  grains.  This  is 
expected  because  slip  is  much  more  difficult  in  the  colony  grains,  particularly  at  lower  strain  amplitudes.  To 
consider  the  orientation  of  the  grains  predicted  to  exhibit  the  estimated  extreme  value  response,  Schmid  factors  were 
identified  for  the  primary  a  HCP  slip  systems  including  the  basal,  prismatic,  pyramidal  <a>,  and  pyramidal  <a+c> 
slip  systems.  In  Figure  7,  the  Schmid  factor  for  each  type  of  slip  is  compared  to  the  extreme  value  volume  averaged 
FS  FIP  for  each  simulation.  In  most  cases,  the  grains  with  the  maximum  grain  averaged  FIP  are  oriented  for  easy 
basal  slip.  Additionally,  prismatic  slip  is  also  observed  but  much  less  often.  When  the  magnitude  of  the  Schmid 
factor  for  prismatic  slip  is  high  (i.e.,  close  to  0.5),  the  magnitude  of  the  Schmid  factor  for  pyramidal  <a>  and 
pyramidal  <a+c>  also  tend  to  be  fairly  high. 

The  extreme  value  marked  correlation  functions  constructed  from  the  simulated  SVEs  of  Ti-6A1-4V,  also  indicate 
that  primary  a  grains  oriented  with  a  Schmid  factor  for  basal  slip  between  0.45  and  0.5  (easy  slip)  are  much  more 
likely  to  be  found  at  the  location  of  extreme  value  response  than  in  the  overall  micro  structure  (see  Figure  8).  This 
agrees  with  what  was  found  by  Bridier  et  al  [15,  16],  who  investigated  fatigue  crack  formation  in  Ti-6A1-4V.  When 
considering  multiple  correlated  micro  structure  attributes,  it  was  much  more  probable  to  find  a  primary  a  grain 
oriented  for  easy  basal  slip  or  prismatic  slip  near  hard  oriented  a+P  colonies  with  high  Schmid  factors  for  BOR 
modified  bee  slip  or  primary  a  grains  with  high  Schmid  factors  for  pyramidal  <a+c>  slip  at  the  locations  of  the 
extreme  value  FS  FIP  than  in  the  overall  micro  structure  (see  Figure  9).  The  latter  previously  mentioned  correlated 
orientation  relationships  were  also  identified  by  Gilbert  and  Piehler  [17]  and  Bantounas  et  al  [42]  in  Ti-6A1-4V. 
Various  work  in  a+P  Ti  alloys  have  demonstrated  that  these  types  of  hard-soft  grain  combinations  can  dramatically 
increase  the  local  stress  state  which  in  turn  promotes  increased  plasticity  in  these  regions  [26,  43-47].  Thus,  similar 
to  that  of  IN  1 00,  these  high  contrast  orientation  relationships  between  softer  basal  or  prismatic  grains  oriented  for 
easy  slip  near  harder  grains  oriented  for  pyramidal  or  bee  BOR  modified  slip  are  predicted  to  significantly  increase 
the  local  driving  forces  for  fatigue  damage  formation. 

5.  Conclusions 

This  work  has  demonstrated  the  use  of  a  new  microstructure-sensitive  extreme  value  probabilistic  framework  that 
links  the  extreme  value  driving  forces  for  fatigue  crack  formation  in  a  P/M  IN  100  and  Ti-6A1-4V  to  the  important 
correlated  microstructure  attributes  that  most  influence  that  extreme  value  response.  Specifically,  we  conclude  that: 

•  Multiple  statistical  volume  elements  (SVEs)  can  be  constructed  and  simulated  using  appropriate  constitutive 
models  to  estimate  extreme  value  distributions  of  key  response  parameters. 

•  Cyclic  plastic  strain  based  fatigue  indicator  parameters  calculated  over  nonlocal  averaging  volumes  on  the  scale 
of  the  microstructure  attributes  important  to  the  mechanisms  of  fatigue  crack  formation  can  be  used  to  identify 
key  microstructure  attributes  important  to  the  extreme  value  response. 

•  In  agreement  with  experiments,  the  simulations  for  IN  1 00  predict  that  cube  slip  may  play  an  important  role  in 
fatigue  damage  formation  particularly  when  there  are  multiple  grains  oriented  for  cube  slip  clustered  in  the  same 
region  or  clusters  of  grains  oriented  for  cube  slip  surrounded  by  other  grains  oriented  favorable  for  octahedral  slip 
[14].  Thus,  as  fatigue  crack  formation  is  dominated  by  the  presence  of  non-metallic  inclusions  or  voids  for  this 
particular  material  system,  we  expect  that  given  two  inclusions  of  similar  character  that  fatigue  cracks  will  form 
preferentially  near  the  inclusion  that  is  surrounded  by  grains  or  grain  clusters  unfavorably  oriented  for  octahedral 
slip. 

•  In  agreement  with  experiments,  the  simulations  for  Ti-6A1-4V  predict  that  primary  a  grains  oriented  favorably 
for  basal  slip  (with  a  Schmid  factor  for  basal  slip  between  0.45  and  0.5)  have  a  much  higher  probability  of  being 
associated  with  the  regions  of  extreme  value  behavior  than  in  the  overall  micro  structure  [15,  16].  Additionally, 
hard-soft  grain  interactions  between  primary  a  grains  oriented  favorably  for  basal  slip  and  primary  a  grains 
oriented  favorably  for  pyramidal  <a+c>  or  a+P  colonies  oriented  favorably  for  bee  BOR  modified  slip  are 
important  relative  to  increased  driving  forces  for  fatigue  crack  formation  [17,  42]. 
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